Performance engineering of semiconductor spin qubit systems 
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The performance of a quantum computation system is investigated, with qubits represented by 
magnetic impurities in coupled quantum dots filled with two electrons. Magnetic impurities are 
electrically manipulated by electrons. The dominant noise source is the electron mediated indirect 
coupling between magnetic impurities and host spin bath. As a result of the electron mediated 
coupling, both noise properties and the time needed for elementary gate operations, depend on 
controllable system parameters, such as size and geometry of the quantum dot, and external electric 
and magnetic fields. We find that the maximum number of quantum operations per coherence 
time for magnetic impurities increases as electron spin singlet triplet energy gap decreases. The 
advantage of magnetic impurities over electrons for weak coupling and large magnetic fields will be 
illustrated. 



PACS numbers: 75.50.-y,75.50.Pp,85.75.-d 

Semiconductor nanostructures [T] are a promising host 
material for the realization of quantum computers [5] be- 
cause of the well developed production technologies and 
because of the potential for scaling to multi qubit sys- 
tems. The greatest bottleneck for realizing semiconduc- 
tor quantum devices is the high degree of noise present 
in the host material. 

Recently substantial progress has been made in the ex- 
perimental demonstration of electron spin coded qubits 
in lateral double QDs [3J, first theoretically proposed by 
Loss and DiVincenzo [I]. This system exhibits promis- 
ing properties for quantum information processing, as it 
can be controlled electrically like charge qubits, however 
offers the longer decoherence times of spin qubits. Nev- 
ertheless, dephasing in spin coded qubits due to nuclear 
hyperfine interaction still presents the major hurdle for 
realizing a scalable quantum computer [5] . 

Another seminal idea of realizing quantum computa- 
tion in semiconductors is by the use of a hybrid qubit, 
first proposed by Kane in a 31 P doped silicon host 
[HOIS]. In silicon, 31 P is a positive dopant that is trans- 
formed into a positively charged nucleus and a loosely 
bound valence electron. Whereas, the qubit is repre- 
sented by the nuclear spin of 31 P, electric gate operations 
on the electron are used for manipulation and diagnos- 
tics of the qubit. The technological implementation of 
this system is more challenging, mostly as a result of the 
significantly smaller dimensions of the qubits. 

We suggest and analyze here a qubit that is a mix- 
ture of the spin coded [4] and of the hybrid qubits 18] 
discussed above. The proposed device consists of a QD 
containing electrons and a neutral dopant acting as a 
magnetic impurity (MI). Similar to the original hybrid 
qubit, the MI represents the qubit, and the electrons 
are used for manipulating the qubit. The use of Mis 
alone for spin-based quantum devices and its noise spec- 
troscopy has been investigated recently [TUl H2] . Here, 
we confine our analysis to Mis with zero nuclear spin, 



[e.g., 56 Fe/ 160 Gd doped in II-VI/IV-VI materials], to ex- 
clude coupling between spins of the electron and nucleus 
of the MI. As shown in Fig[TJ the MI is localized in space 
and its direct interaction with the host spin bath is neg- 
ligible. Therefore, the decoherence of the qubit is deter- 
mined by the electron mediated coupling of the MI to the 
spin bath. 

Our configuration has two advantages over the original 
qubit systems. First, it allows greater design flexibility, in 
particular with regard to the size of the confining poten- 
tial. This will facilitate the actual technological realiza- 
tion of a hybrid qubit. Second, Mi-electron coupling and 
electron mediated Mi-spin bath coupling show particular 
dependence on the electron spin singlet triplet energy gap 
A e , that allows forming stable qubit over a range of QD 
confining potentials, and the external electric and mag- 
netic fields. This system corresponds to the electron spin 
coded qubit system investigated in Ref. [4] and makes a 
comparison possible. The performance of our quantum 
computing setup is measured by the maximum number of 
operations Ni — Ti/Ti with i = e, m for electrons and Ml 
qubits, respectively. Here, r is the coherence time and 
T is the time required for an elementary gate operation. 
In this work we illustrate that the quantum performance 
of the spin coded qubit based on Mis increases by the 
external electric and magnetic fields. This is in contrast 
to the electron spin coded qubit that its performance is 
suppressed rapidly by the inter-dot coupling. 

We represent the QD system by the Hamiltonian H = 
H c + H m + H n + Hi n , which contains electrons (H e ), Mis 
(iJ m ), host semiconductor nuclei (H n ), and their inter- 
actions (H- ln ). In our model we consider e-MI and e-n 
interaction, i.e. H- m = H crn + H cn . The direct interac- 
tion between MI electrons and semiconductor host nuclei 
are neglected, as d- and f-electrons are highly localized in 
space. Further, our analysis is confined to Mis with zero 
nuclear magnetic moment, which do not interact with 
the MI electrons. As a result, dipole-dipole interaction 
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FIG. 1: Schematics of two spin qubit unit in a DQD consid- 
ered here. Each spin qubit is represented by a MI. Interaction 
between qubits is mediated by the two electrons; the electron 
MI coupling coefficient is denoted by J em . As interaction of 
Mis with nuclear spin bath is negligible, the dominant noise 
source is electron mediated coupling between Mis and nuclear 
spin bath; the electron nuclear spin bath coupling coefficient 
is given by A. The interaction strength between the Mis and 
between MI and nuclear spin bath depends on the singlet- 
triplet electron energy gap A e , which can be controlled by 
external electric gate voltage (V 9 ), magnetic field (B), and 
shape/size of the DQD. This allows active noise engineering 
and optimization of the two qubit system. 



between the nuclear spins of MI and host nuclei has no 
influence on our system. 

Electrons confined in quasi-two-dimensional quantum 
dots in a uniform perpendicular magnetic field can be 
described by the effective mass Hamiltonian 



(ifaa) are calculated by discretizing T + Z in real space, 
and diagonalizing the resulting matrix. By using the cre- 
ation (annihilation) operators c' arT (c aa ) for an electron 
in a non-interacting single-particle (SP) state \a,a), the 
Hamiltonian of an interacting system in second quanti- 
zation can be written as 

H W , c t c 
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where the first term represents the single 
particle Hamiltonian and V aa ^ a i tVa = 
J drj d? V l a {f^ip^,{?) 7 ^^^>{?)ip va {T f ), is the 

two-body Coulomb matrix element. 

The Hamiltonian for the Mis accounts for MI-MI direct 
exchange interaction and Ml-Zeeman coupling, 

M 

H m = J ii'Mj ■ Mf +J29m»hM zj B, (3) 
i,j'=i j 

where Jjy is the direct MI-MI antiferromagnetic (AFM) 
coupling, g m is the MI g-factor, and M Z j is the z- 
component of the MI spin operator. 

The nuclear-nuclear direct dipole interaction in the 
host semiconductor is neglected. The nuclear Hamilto- 
nian is given by the Zeeman coupling term, 



(4) 
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where g n is the nuclear g-factor and I z i is the z- 
componcnt of nuclear spin operator. 

The e-MI exchange interaction is modeled by 



H em = - J cm Y & ■ M]$( r t - R-i)) 



(5) 
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where T = l/(2m*) ((h/i)V + (e/c)A(rf) +V{r) is the 

single electron Hamiltonian in an external magnetic field 
B = Bz, perpendicular to the plane of the confining po- 
tential. Here (f) = (x, y, z) describes the electron posi- 
tion, V(r) denotes the quantum dots confining potential, 
and A(f) = (1/2)5 x r is the vector potential. Further, 
to* is the conduction-electron effective mass, — e is the 
electron charge, and e is the host semiconductor dielec- 
tric constant. Finally, Zi = (l/2)g c ^S Z iB determines 
the Zeeman spin splitting, g c is the electron g-factor in 
host semiconductor, /ib refers to the Bohr magneton, and 
S z i represents the z-Pauli matrix of electron i. 

The single particle eigenvalues (e arT ) and eigenvectors 



with J em the exchange coupling between electron spin 
Si at Vi and impurity spin Mj located at the position 
Rj [10] . In second quantization it can be written as 



Hem = -/2/2 



a/3 I 



(6) 



where J Q ^(Rj) = J em iy9* (Rj)^(Rj). Similarly, we de- 
scribe the electron - nuclear spin bath hyperfine interac- 
tion by 



H ea = ^2A l S i -Ii5(r i -Ri) 



(7) 



i,l 



with Ai = (16n/ 3) HbfJ>i/Ii the isotropic (Fermi contact) 
part of the electron-nucleus hyperfine interaction [BJ 1131 
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[T3] . Here hi and R; are magnetic moment, and position 
of the Ith nucleus and sum goes over all nucleus in the 
lattice. In second quantization it can be written as 



a/3 I 



(8) 



where A a/3 (Ri) = Anp* a (Ki)<p (Ri). Finally 

Hint 2 ^ ; ^ ; Qaa' ^Vcr' C Q ( j C ct'cr / ) (9) 

ftd' crcr' 

where Qaa' = Ej J<xa'(Rj)Mj - E„ A aa >(Ri)Ii. 

From the total Hamiltonian an effective Hamilto- 
nian is obtained by tracing over the degrees of freedom of 
the electron wavefunction and by taking the interaction 
term H ln into account to second order of perturbation 
theory, which yields 



H, 



eft 
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(10) 



Here we limit our calculation to a two electron and two 
MI system in a DQD. The two-electron wavefunction is 
confined to the Hilbert sub-space constructed from the 
bonding and anti-bonding (HOMO, LUMO) one-electron 
orbitals of the DQD, ip±. Below the magnetic field cor- 
responding to spin singlet-triplet transition, this results 
in six basis functions of two-electron, a spin singlet (So) 
ground state Vfg that can be expressed as superposition 
of ip + (ri)ip + (r 2 )\So) and (p_(ri)ip-(r2)\So) with binding 
energy E g , and five excited states \l/ x with energy E x , 
consisting of three degenerate first excited triplet states, 
and two higher excited singlet states [TT] . 

Calculating the matrix elements in Eq. [10] yields the 
effective Hamiltonian 



H c ff — H n 



(11) 



where H mm = Y)j ji A jy-Mj ■ My is the electron medi- 
ated (RKKY-type [M]) interaction between the Mis, and 
H mn = Ej i Ajlll ' Adj is the electron mediated interac- 
tion between Mis and nuclear spin bath; the electron 
mediated interaction between host nuclear spins H na is 
neglected. Here, A,y = -j 2 Jl n U(R :j ,R jl )/(2A e ) and 
A jt i = ^ 2 AJ em U(R v Ri)/l\ e . Further, 7 = 0+- a_, 
where a+ , and a_ are the coefficients of the two-electrons 
ground state that is expressed as linear combination of 
bonding- antibonding in two level model ^(fj^fa) = 
[a + Lp + (ri)Lp + (r 2 ) + a-ip-(fi)(p-(f 2 )] \S Q ), A e is the two 
electron singlet-triplet splitting, A = 1/iEfci^Zs and 
U(Rj,R x ) = v+{R 3 )^Rj)v+{R\)V-{R\) with A - 

In the following we use the effective Hamiltonian [TTJ to 
calculate the decoherence time of the DQD with two Mis 
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FIG. 2: iVi, the maximum number of elementary gate op- 
erations within the coherence-time normalized by factor 2-kA 
(A = A e « A m ) in meV in a FeiZnSe DQD (2MI), com- 
pared with ZnSe DQD filled with two-electrons (2e) versus 
parabolic strength of confining potential of each dot, u>o, for 
various gate voltages, V 3 — 110 (circles), V 3 — 154 (squares), 
V g = 198 (triangles), and V g = 242 (diamonds) in meV. Ni as 
a function of external magnetic field for ujq — 27.5 meV and 
V g = 110 meV is shown in the inset. 



and two electrons, where each MI represents a spin qubit. 
The result is compared to a two qubit system realized by 
two electrons in a DQD. The calculation is performed by 
using the quasi-static bath approximation |15H17j , where 
the host nuclear spins are approximated by a random 
magnetic field B n with a Gaussian distribution. In this 
limit the two-electron and two-MI nuclear bath Hamilto- 
nian are given by 



#kn — ^2 5kMb-B n • Ki 



(12) 



i=i 



where k = e, m and K = S, M for electrons and MI, 
respectively. The coherence time is obtained by solving 
the equation of motion for K\ and K 2 with initial state 
|tl) and by averaging over the Gaussian magnetic field 
distribution. From that we obtain (B c ) = (B m ) = 0, 
(Bt) = l/GfewO'E,/^ + l)if|^(R/)| 4 , and (Bl) = 
l/(.9mMb) 2 E/ + 1)^ 2 |<P+(RzV-(Rz)| 2 - From there 
an effective Zeeman splitting Ak = (2(B 2 ) /3) 1 / 2 is cal- 
culated, hence Tk = ^/(gk/^bAk). Assuming = 1/2 
we find the spin relaxation time t c = h/(2A e ) and 
r m = Ch/(2A m ). Here A e = (Ez ^+(R;)| 4 ) 1/2 , 
A m = (Ei^l^+CR-O^-CROI 2 ) 172 , and C = 
A c /[7 2 J Gm A(fli,-R 2 )] is the RKKY correction to the 
MI coherence time, stems from the Ml-nuclear-spin in- 
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teraction mediated by electrons. Here A(R 1 ,R 2 ) — 

\ l P+(^j) l P-(^'j)\ describes the spatial dependence 
of the e-MI exchange interaction, a parameter that de- 
pends on the electron envelop wavefunction at the MI 
positions Rj. Note that in the limit of zero inter-dot 
tunneling, A e = A m . The ratio of MI and electron co- 
herence times can be calculated as 



I™ 

To 



A, 



U em A(R 1 ,R 2 )A r , 
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(13) 



The performance of a quantum computing setup is 
given by the maximum number of operations iVj = Tj/Tj 
with i = e, m for electrons and MI qubits, respectively. 
Here, r is the coherence time and T is the time required 
for the elementary gate operations. Our system is com- 
pared to the original proposal in Ref . [4] , where an XOR 
gate control in a two-electron DQD is analyzed. The el- 
ementary gate operations needed for the XOR gate are: 
(i) a correlated spin swap from — >|-|t)> where the 
first and second position refers to the left and right dot, 
respectively, and (ii) single qubit operations with an ex- 
ternal pulsed magnetic field. The second operation has to 
be done within the time of one spin swap. Therefore, the 
time for one XOR gate is determined by the correlated 
spin swap time T. 

The MI and electron spin swap times are given by 
A" 1 , and A^ 1 modulus nh, assuming that the gate 
voltage and therewith exchange coupling is controlled 
by a square pulse [1]. Here, A m = Ajj> is the cou- 
pling coefficient between the two Mis. As a result, 
N c = A c /(2irA e ), N m = [J em Z7/(2A)]/(27rA m ) where 
U = p + ( J 3i)p_(£ 1 )¥> + $ a )¥J_( J £ 2 ) and R u R 2 are MI 
coordinates. The ratio of the maximum number of el- 
ementary operations (per coherence time) is given by 
N m /N e = (r m /r e ) (A m /A c ), which finally gives 



^em U A e 

A r 2AA„ 



(14) 



In Eq. (14 1, N m /N e cx 1/A C . Unlike the other param- 
eters in Eq. (14 1, A c decays to zero very rapidly by in- 



creasing the external magnetic field and inter-dot energy 
barrier that lowers the inter-dot coupling. Therefore one 
expects to observe decay in performance of the electron 
spin coded qubit due to variations in B and V g . Un- 
like the electrons, Mis show a robust increase in their 
quantum operation performance. To gain the optimum 
performance of Mis over electrons we employ a numerical 
calculation based on exact diagonalization of Eq. ([T]) from 
which the input parameters for Eq.( 14) are obtained |llj . 
The DQD is chosen to be double Gaussian along the axis 
and parabolic in the perpendicular direction and each MI 
is centered at one of the QDs. We perform our calcula- 
tion for three different materials, Gd:PbTe, Fe:CdSe, and 
Fe:ZnSe where they show T m /r c 300, 20, 2 respectively 
at V g = 154 meV and B = 0. The values for J cm are 



adopted from Ref. [IS]. In this range of parameters the 
coupling between Mis and nuclear spins (Ajj) is opti- 
mized to be weak to maximize r m /r e . However with de- 
creasing Aj i, the coupling between two Mis (A m ) lowers 
and as a result the time required for fundamental gate 
operations becomes longer. Thus, the gain in coherence 
time (r m /T c ) is offset by a loss in gate operation times 
(A /A m ). To maximize N m /N e we search for a range 
of parameters that allows simultaneous maximization of 
T m /r e and A m /A c . 

In Fig. [2]we show the results for Fe:ZnSe. N m and N e 
normalized to 2irA e are plotted as a function of parabolic 
confining strength ljq for various values of the gate volt- 
age V g , and the external magnetic field B (inset) . Within 
numerical parameters considered in this calculation we 
found A e w A m . We observe that N m /N e increases 
with increasing B, V g (inter-dot energy barrier) and coq 
(tighter confinement). Within the parameter range con- 
sidered here, a maximum performance increase of about 
three orders of magnitude can be achieved over the elec- 
tron spin coded qubit. It is important to mention that 
even higher increases of N m /N e might be achievable. The 
maximum B and V g values used for our optimization had 
to be limited to the range of validity of our two level 
model. For increasing values of B and V g the system 
approaches the singlet triplet transition point at which 
A c — > and N m /N c cx 1/A — >• oo. In this limit our 
approximation based on the two level model fails, and 
a more exact analysis becomes necessary. This will be 
studied in more detail in a follow-up work. Further, the 
smallest realizable quantum dot size is around 5 nm cor- 
responding to the maximum luq in Fig. [2] However, an 
extrapolation of the numerical results to the atomic scale 
indicates performances gains N m /N e of more than five 
orders of magnitude. In this limit, our system becomes 
comparable to Kane's proposal [5J, which demonstrates 
its favorable performance properties. The advantage of 
our system is that a compromise can be found between 
optimizing performance and accommodating technologi- 
cal limitations. 

In conclusion the qubit system investigated here opens 
the possibility for noise and performance optimization. 
We have found that a combination of active and passive 
optimization is necessary to obtain appreciable improve- 
ments; our analysis of the hybrid Ml/electron qubit pre- 
dicts a performance gain of at least three orders of magni- 
tude over electron spin coded qubits in the limit of small 
inter-dot coupling. 
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